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Abstract. Liu procedure is applied to a first order weakly nonlocal special 
relativistic fluid. It is shown, that a reasonable relativistic theory is and ex- 
tended one, where the basic state space contains the momentum density. This 
property follows from the structure of the energy-momentum balance and the 
Second Law of thermodynamics. Moreover, the entropy depends on the en- 
ergy density and the momentum density on a given specific way, indicating 
that the local rest frame energy density cannot be interpreted as the inter- 
nal energy, the local rest frame momentum density should be considered, too. 
The corresponding constitutive relations for the stress and the energy flux are 
derived. 



1. Introduction 

Nonrelativistic nonequilibrium thermodynamics separates the dissipative and 
nondissipative parts of the evolution of physical quantities. This separation is 
based on the construction of the internal energy balance [2 [51 [3] . According to 
the classical interpretation the internal energy is the difference of the total energy 
and the known specific energy types. The entropy function depends directly on 
the internal energy. The internal energy is distributed unbiased among the molec- 
ular degrees of freedom. The process how the other energy types are converted to 
internal energy is called dissipation. This approach is common in every theories 
of nonequilibrium thermodynamics including classical irreversible thermodynamics, 
where the hypothesis of local equilibrium applies. However, there is nothing simi- 
lar in relativistic irreversible thermodynamics. Furthermore, practically there is no 
relativistic irreversible thermodynamics at all, because the local equilibrium theory 
is plagued by serious controversies, therefore only the extended theories, theories 
beyond local equilibrium, are considered as viable. The reason is that the classical 
theory of Eckart for relativistic fluids \A\ is simple and elegant, but produces generic 
instabilities [5]. The more developed extended theories incorporate the theory of 
Eckart, therefore inherit (but more or less suppress) the instabilities [6l[7l[8j. 

In this paper we investigate the possibility of local equilibrium in relativistic 
hydrodynamics by methods of continuum thermodynamics. At the next section 
the balances of energy-momentum and entropy are introduced. In the third section 
we calculate the dissipation inequality for local equilibrium (first order) relativistic 
hydrodynamics by Liu procedure. The need of second order (extended, or weakly 
nonlocal) theories is indicated by the emergent structure. A new concept of rel- 
ativistic internal energy follows. Based on these results we give the constitutive 
equations of the simplest extended theory by the heuristic arguments of irreversible 
thermodynamics in the fourth section. 
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2. Basic balances of relativistic fluids 

For the metric (Lorentz form) we use the g^"^ = diag{—l, 1, 1, 1) convention and 
the units are introduced that the speed of hght c — 1. Therefore for a four- velocity 
we have Uau" = —1. A"^ = g"^ + u^up denotes the w-orthogonal projection. 
First we give the basic balances of energy-momentum and entropy. 

The energy-momentum density tensor is given with the help of the rest-frame 
quantities 

(1) T"^ = eu"M^-f M"/-|-u''g" + P"'', 

where e = UaUf^T"-^ is the density of the energy, — —Ua^^jT"'' is the energy 
fluxoT heat flow q" — —UfjA'^^T''^ is the momentum density and P"^ = A^A^T'^^ 
is the pressure tensor. The momentum density, the energy flux and the pres- 
sure are spacelike in the comoving frame, therefore Mq?" — and Uaq" = and 
UaP"^ — UaP^" = 0^. The energy-momentum tensor is symmetric, because we 
assume that the internal spin of the material is zero. In this case the energy flux and 
the momentum density are equal. Let us emphasize that the form ([1]) of the sym- 
metric energy-momentum tensor is completely general for one-component fluids, 
nevertheless it is expressed by the local rest frame quantities. 

Now the conservation of energy-momentum dpT"^ = is expanded to 

(2) dfsT'^'^ = eu" + eu'^dpu'^ + eif- M"9/3g'^ + /9/3m" + + g"9/3u^ + 9^3^"'', 

where e = = u'^daG denotes the derivative of e by the proper time t. Its 
timelike part in a local rest frame gives the balance of the energy 

(3) - UaS^jT"^ = e -f ec)„u" + d^q"" + q^^ii^ + P^'^dpUa = 0. 

The spacelike part in the local rest frame describes the balance of the momentum 

(4) A%dpT^P = eit" + q^dpuf^ + q^dpu" + A%q< + A'^^dpP^" = 0". 

The entropy density and flux can also be combined into a four-vector, using local 
rest frame quantities: 

(5) 5*" = 81*" + J", 

where s = —UaS" is the entropy density and J" = S*" — = A^^S^ is the 
entropy flux. The entropy flux is u-spacelike, therefore UaJ" = 0. Now the Second 
Law of thermodynamics is expressed by the following inequality 

(6) SaS*" = S + s9aW" +aaJ" > 0. 

3. Thermodynamics 

The thermodynamical background in relativistic theories is usually based on 
analogies with nonrelativistic thermostatics. However, nonequilibrium thermody- 
namics developed beyond the simple 'let us substitute everything into the entropy 
balance and see what happens' theory since Eckart. It is important to check the 
dynamic consistency of the Second Law, considering the evolution equations as con- 
straints for the entropy balance. This method of nonequilibrium thermodynamics is 
constructive, gives important information for new theories and reveals some deeper 
interrelations. Here we exploit the Second Law by Liu's procedure [9] introduc- 
ing a first order weakly nonlocal state space in all basic variables, thus restricting 
ourselves to a local equilibrium theory. One can find a general treatment of nonrel- 
ativistic classical and extended irreversible thermodynamics from this point of view 
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in [lOj . Our aim here is to investigate the relativistic fluids with similar methods 
to get the relativistic equivalent of the classical Fourier-Navier-Stokes system of 
equations for one component fluids. 

Our most important assumption regarding relativistic thermodynamics is that 
the constitutive equations are local rest frame expressions. From a physical point 
of view it is natural, because material interactions are local. 

The basic state space of the theory is spanned by the energy density e and by 
the velocity field u" . The constitutive state space is spanned by the basic state vari- 
ables and their first derivatives, therefore it is first order weakly nonlocal. Hence 
the constitutive functions depend on the variable set C — {e, Ua, daC, daUp). The 
constitutive functions are the energy flux/momentum density q" , the pressure P"'', 
the entropy density s and the entropy flux J". The derivatives of the constitutive 
functions are denoted by the serial number of the corresponding variable in the 
constitutive space, e.g. g^g^^^ = 93s. With this notation we can distinguish eas- 
ily the derivatives by the constitutive and spacetime variables. A nonequilibrium 
thermodynamic theory is considered to be solved if all other constitutive quantities 
are expressed by the entropy density and its derivatives. 

According to the procedure of Liu the balance of energy-momentum (|2| is a 
constraint to the entropy balance ^ with the Lagrange-Farkas multiplier A^ 

(7) ^a^" - A^dpT"'^ > 0. 

Let us remember, that here the spacelike components of the four quantities and 
the entropy density are the constitutive quantities depending on the introduced 
constitutive variables C. Therefore, in the above inequality we can develop the 
derivatives of the composite functions. The coefficients of the derivatives that are 
not in the constitutive space must be zero, therefore we get the following Liu- 
equations: 

(8) da^pe : (dsS'^f -A^idsT^^ = 0"", 

(9) da^pu^ : {diS'^ f'' - Af.idiT^'^f'' = 0'^'^^. 

The simple structure of the Liu equations suggests the assumption that the 
Lagrange multiplier is a local function, it does not depend on the derivatives of the 
basic state variables 

(10) A^ = A-f{n,e). 
Then a general solution of (I!])-® is 

(11) S" - A^T'>°' - A" = 0", 

where A" = A" (n, e) is an arbitrary local function. 

Let us introduce the splitting of the vector multiplier and the four-vector A" to 
spacelike and timelike parts in the local rest frame as 

A" = -Au" + r, 
where for the spacelike components UqP = Uao" = 0. Now equation (jlip gives 



(12) m"(s - Ae - l^q'' -A) + (J" - Ag" - Z^P'^" - a") = 0". 
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Here both the timehke and spacehke parts are zero, resulting in 

(13) s = Ae + l-yq'f + A, 

(14) J" = Ag" + Z^P'^" + a". 

After the identification of the Liu equations we can get the dissipation inequality 

as 

dae [{dis)u" + di.r - Au" - ASig" - ;^9iP^" - Z^^ig^u"] + 
d^ui3 [{s - Ae - l^q^) A"^ + {d2sfu" + (^2 J")''- 

(15) if^eu" - If^q" - k{d2q^T - A^{d2P''°'f - A^u"(92g'')^] > 0. 

Here we exploited that the partial differentiation by e can be exchanged with 
the multiplication by the four velocity u". 

In the dissipation inequality one should consider the solution of the Liu-equations. 
Substituting ((T31) and into ((15]) we get 

9„e [(9is - A - l^diq'')u°' + q^diA + P^°'dil^ + Si a"] + 

d^up [AA"^ + q°'{d2h-Y + P'<''{d2l-yf + {d2a''f+ 

(16) {{d2sf - l^{d2q''f ~ fe ~ Aq^) - l^^q'^ - AP'^^] > 0. 

Here the following identities were applied simplifying the last term (92 = duf^) 

u^duf,q^ = du^{u-^q'^) - q'^duf.u^ = ~q'' = ~q'^ , 
u^du.P'"' = du, {u^P^^) - P^'^du.u^ = -P-'^A/ = -P^". 

Observing the first term in the last form of the dissipation inequality one can 
eliminate the direct velocity dependence of the entropy function recognizing that 
the entropy may depend on the energy flux in the following form 

(17) s(e, u^jdaC, dau'^) = s{e, q''{e, u",dae, daU^))- 

Therefore the entropy is local, independent of the derivatives of the basic state 
space variables and the velocity field, but depends on the energy flux. The energy 
flux can depend also on the derivatives because according to our initial assumptions 
it is considered as a constitutive function. Then the Lagrange-Farkas multipliers 
are determined by the entropy derivatives 

(18) deS = A, dqcS = la- 

We introduce a temperature T as 

(19) dJ = A=^. 

We may recognize that a full thermostatic compatibility requires that in (|13p : 
A := where p is the pressure. These consequences are completely analogous 
to the results of the nonrelativistic nonequilibrium thermodynamic theory, where 
thermostatics arises from the structure of the balance form evolution equations as 
constraints for the Second Law. 

Finally we assume that entropy fiux is classical and the additional term a" [11] 
is zero 

(20) a" = 0". 
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Then the dissipation inequahty reduces to the following simple form 

(21) - i (P"/' + T/^" - pA'^") d^up - P^-^dJ^ ^ ( Y + ^'") " 

As we do not want an acceleration dependent entropy production, we must re- 
quire that the last term vanishes. According to (fT8|) and (fT9| 

(22) edqJ + q°'deS = Q. 
The general solution of can be given as 

(23) s = S(e2-g"g„)+B, 

where i?=const. The entropy must depend on the energy density e and the mo- 
mentum density g" on a very particular but simple way. As a consequence of this 
functional form of the entropy function the Gibbs relation can be given with the 
help of the entropy derivatives pS]) as 

(24) de - —dqc = Tds. 

e 

We may require first order homogeneity of the entropy in psp , without restricting 
the generality. That can we get introducing E = \/\e^ — qaq" as a variable of the 
entropy density. In this way the entropy is a first order homogeneous functions 
both of the energy density e and the momentum density q". With this property it 
is unique. 

The corresponding potential relation can be constructed according to the first 
order homogeneity of the physical quantities as 

(25) = Ts-p. 

e 

The previous thermostatic relations require the interpretation of E as internal 
energy. On the other hand let us recognize that E is the absolute value of the 
energy vector 



(26) E = \\E^ = II - upT^'^W = \\eu" + q'^\\ = VW^Q^l 

However, one should pay attention that the l/T introduced in (fT9)) is not the 
derivative of the entropy function according to E. 
Finally the entropy flux from (fT4|) and ([20|) 

(27) J" - ^g" - ^P^". 
The final form of the dissipation inequality is 

(28) g"9a ^ - ^ (^P"^ + pA"^^ do^up - P^-'da. 1^ > 0. 

The last term in this expression with a derivative of one of the constitutive 
quantities indicates that we cannot give proper thermodynamic fluxes and forces, 
as a solution of the inequality. An other problem appeared already with (P^ . 
because /„, as the spacelike part of the Lagrange multiplier in a local rest frame, 
was assumed independent of the derivatives of e and . That is the Fourier heat 
conduction was excluded as a possible constitutive function. Both problems indicate 
that a complete theory may exist only either in an enlarged constitutive space or in 
an extended basic state space. One of the possibilities is to introduce higher order 
derivatives of the basic state space into the constitutive state space and construct 
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a second order weakly nonlocal theory. The other possibility is to enlarge the basic 
state space and construct an extended theory. In both cases the key that may lead 
beyond the traditional Miiller-Israel-Stewart theory is the new internal energy E. 



4. Extended irreversible thermodynamics of relativistic fluids 

Motivated by the results of the previous section we calculate the entropy produc- 
tion by a direct substitution of the balance of the energy into the entropy balance. 
We are to construct an extended theory, introducing as an independent variable, 
but exploiting the fact that the entropy depends both on the energy and momentum 
densities by the specific way derived above. 

The entropy flux is assumed to have the very classical form 

(29) J" = ig". 

Substituting the energy balance ^ into the entropy balance equation we arrive 
at the following entropy production formula: 

-^(ea„u" + 9„g" + + P"^5^u„) + |^g„ + sS^u" + (^^g"^ = 

(30) {P"^" i-e + sT)A^'') d^up + g" - ^ - ^) > 

In isotropic continua the above entropy production results in constitutive func- 
tions assuming a linear relationship between the thermodynamic fluxes and forces. 
The thermodynamic fluxes are the viscous stress 11"'^ = (P"'^ — (sT — e)A"''), and 
the energy flux q". We get 

(31) n"^ = P"^ - A"'' (^p - = -2r/(A"''A''^5^u^)"° - r/^5^w^A"'^, 

(32) g" = -X^A''^ (d-,T + Tu'' + ^ 

i ^ \ e 

where denotes symmetric traceless part of the corresponding second order tensor 
(A^^y^ = ^{A^^ + A^^) — ^A'^S^^ and we have introduced the scalar thermostatic 
pressure according to (^5]) (therefore p ^ P"/3). (Plj) and 1(5^ are the relativistic 
generalizations of the Newtonian viscous stress function and the Fourier law of 
heat conduction. The shear and bulk viscosity coefficients, rj and rj^ and the heat 
conduction coefficient A are nonnegative, according to the inequality of the entropy 
production ([30]) . 

The equations ([3]), (jH) are the evolution equations of a relativistic heat conduct- 
ing ideal fluid, together with the constitutive functions ([5T]) and As special 
cases we can get the relativistic Navier-Stokes equation substituting (|3T|l into (|4]) 
and assuming q" ~ 0, or the relativistic heat conduction equation substituting (|32p 
into Q assuming that 11"'^ = 0. The heat conduction part results in a special 
extended theory, where only the energy flux appears as independent variable. 
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5. Summary and discussion 

Int the first part of the paper we have investigated the local equilibrium theory of 
special relativistic fluids. We have seen that there may be no such theory that could 
give a complete solution of the entropy inequality provided the following conditions 

(1) local Lagrange- Far kas multipliers, 

(2) local entropy p7|) . 

(3) no additional term in the entropy flux ((20)) . 

The first two assumptions were necessary to get a particular solution of the Liu 
equations and the dissipation inequality. On the other hand they are natural in 
local equilibrium. 

We have concluded that either an extension of the basic state space or an enlarge- 
ment of the constitutive state space can give a complete solution. Our investigations 
indicated a particular dependence of the entropy on the energy and momentum den- 
sities, leading to a distinction of internal and total energy densities of relativistic 
fluids. 

The local rest frame energy density e = UaT^^up is usually interpreted as in- 
ternal energy in thermodynamic theories. However, the symmetry of the energy- 
momentum tensor can hide that energy flux is related to dissipation, but momentum 
density is not. According to the previous investigations the total energy density e 
- the time-timelike part of the energy-momentum tensor - is not a suitable internal 
energy, the entropy density should be a function of the absolute value of the energy 
vector i5" = —upT°'^ - the timelike part of the energy momentum. 

To compare our proposal to the traditional Miiller-Israel- Stewart theory [l2llT3] 
it is instructive to expand the internal energy into series, assuming that > q^qa 

(33) S= v/|e2-(?"(Z„| «e-^ + .... 

2e 

The last, quadratic term in the above expression is what appears in the Miiller- 
Israel-Stewart theory. However, in our case 

- the corresponding relaxation time is fixed r — 1/e, 

- the quadratic term is only the first approximation, 

- only the energy flux was introduced as an independent variable in our ex- 
tended theory, the viscous stress was not necessary. 

The series expansion is instructive comparing to nonrelativistic hydrodynamics. 
There the internal energy is the difference of the total energy and the relative kinetic 
energy. In ([55]) the quadratic expression is what one could consider as a kind of 
energy of the flow, but of course it is not connected to an external observer, it 
considers only the local rest frame momentum density. In a sense our expression 
shows that introducing E as internal energy we declared that the energy of the flow 
(in the local rest frame) does not give a dissipative contribution. 

The extension of the present calculations considering the balance of particle num- 
ber is straightforward. Moreover, one can show that the above system of equations 
gives a stable homogeneous equilibrium in linear stability investigations, contrary 
to the theory of Eckart [M], therefore it can be considered as a minimal viable 
extension of the local equilibrium theory without the complexity of the Miiller- 
Israel-Stewart one. 
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